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Abstract 

We consider the problem of the speed selection mechanism for the one di- 
mensional nonlinear diffusion equation Ut = u xx + f(u). It has been rigorously 
shown by Aronson and Weinberger that for a wide class of functions /, suffi- 
ciently localized initial conditions evolve in time into a monotonic front which 



propagates with speed c* such that 2^J /'(0) < c* < 2-\/sup(/(u)/u). The 



lower value cl = 2^//'(0) is that predicted by the linear marginal stability 

speed selection mechanism. We derive a new lower bound on the the speed of 

the selected front, this bound depends on / and thus enables us to assess the 

extent to which the linear marginal selection mechanism is valid. 
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I. INTRODUCTION 



In several problems arising in biology, population dynamics, pulse propagation in nerves, 
crystal growth, fluid flow and others, it is found that if the system is suddenly made un- 
stable, the subsequent dynamics is characterized by the propagation of fronts. The systems 
for which this phenomenon occurs have received much attention recently, especially related 
to the problem of pattern formation. A small perturbation at a localized point grows to 
eventually cover the whole space. An important problem to be solved has been the deter- 
mination of the speed at which the front of the pattern moves into the undisturbed regions 
of the system and the wavelength of the pattern left behind. (For a recent and extensive 
review of this subject we refer to [U and references therein.) Several authors @-[5|] have 
formulated criteria that provide an answer to these questions. These criteria are heuristic 
extensions to higher order equations of rigorous results and heuristic arguments which have 
been developed for the nonlinear diffusion equation 



where f(u) 6 C 1 [0, 1], f(0) = f(l) = 0. In what follows we assume that / is positive in 
(0,1). In this case u = is the unstable fixed point and u = 1 is a stable fixed point. 
Aronson and Weinberger || have shown that any positive initial condition Uq(x) < 1 for all 
x, which decays exponentially or faster at infinity will evolve into a front propagating with 
speed c*. This asymptotic speed is the lower speed for which equation (1) has a monotonic 
front joining the stable state u — 1 to the unstable state u — 0. Moreover, 



For the special case of the Fisher-Kolmogorov equation f(u) = u — u 3 , f'(0) = 1 and 
sup(f(u)/u) = 1 so that c* = 2. In general 0, for any concave f(u), sup(f(u)/u) = f'(0), 



and c* = 2y/'(0). The value c* = 2 is the value which had been derived by Kolmogorov, 
Petrovsky and Piskunov |J using an heuristic argument (the linear marginal stability mech- 
anism) which is equivalent to the conjecture that the asymptotic speed of the front is that 
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u t = u xx + f(u) 



(1) 




(2) 



for which a perturbation to the front is marginally stable in the frame moving with the 
front speed. Based on the applicability of this argument for the Fisher Kolmogorov equa- 
tion and more generally for concave functions / several authors have developed extensions 
of this argument to higher order equations. These generalizations are purely heuristic, the 
only rigorous results available being those of Aronson and Weinberger. In general however, 
sup (/(«)/«) is not /'(0) and equation (2) gives a bound on the selected speed. It is known 
that for some choices of /, and explicit examples have been given, c* is greater than 2. 
These cases, referred to as those in which a nonlinear marginal stability mechanism oper- 
ates, have been generalized || for higher order equations based on the observation that for 
the nonlinear diffusion equation the selected front is that with the steepest decay to zero. 
The exact point of transition from the linear marginal stability to the nonlinear regime has 
been determined for functions / of the form f(u) = fiu + u n — u 2n ~ l for which an exact 
solution for a monotonic front can be given. It has been shown that for \x smaller than a 
critical value the solution corresponds to a nonlinear marginal stability solution || . To the 
best of our knowledge, the only lower bound on the speed that, for general /, shows that the 
linear speed is not always preferred, has been given recently by Berestycki and Nirenberg 
T0|| . They show that 



c 2 > 2 [ f(u)du (3) 
Jo 

from where it is evident that for sufficiently large / the speed exceeds the marginal value 
C£. The purpose of this work is to give a new bound that enables one to evaluate the regime 
of validity of the linear marginal stability mechanism with increased accuracy. As shown 
by Aronson and Weinberger, the asymptotic speed of the front is the lowest for which there 
is a monotonic travelling wave solution u = q(x — c*t) of equation (1). The selected speed 
satisfies q zz + c*q z + f(q) = 0, lim , u 2 _ > _ 00 = 1, limu^oo = 0, where z = x — c*t. We find 
it convenient to work in phase space, where monotonic fronts obey an equation of an order 
less than the original equation. Since the selected speed corresponds to that of a decreasing 
monotonic front, we may consider the dependence of its derivative dq/dz on q. Calling 



p(q) = —dq/dz, where the minus sign is included so that p is positive, we find that the 
monotonic fronts are solutions of 

P(q)f q -c*p(q) + f(q) = 0, (4a) 

with 

p(0) = 0, p(l)=0, p>0 in (0,1). (4b) 

The bound follows in a simple way from equation (4a). Let g be any positive function in 
(0,1) such that h = —dg/dq > 0. Multiplying equation (4a) by g/p and integrating with 
respect to q we find that 

l\ hp+ ir 9 ) dq = c * Jo gdq (5) 

where the first term is obtained after integration by parts. However since p, h, f, and g are 
positive, we have that for every fixed q 

f(q)9 



hp+ iSHJL>2^fgh 
p v 



hence we obtain our main result, 



c* > 2 thlpMl (6a) 
Jo gdq 

where 

g > and h = -g' > in (0, 1). (6b) 

That this result yields a better bound than that given by equation (3) can be seen by 
choosing g so that g h = f, and g{l) = [12| . 



Next we illustrate the use of this bound by applying it to two explicit forms of /. Since 
here we wish only to illustrate the use of this bound, we take three simple trial functions. 
As a first trial function choose g so that f = h = —g[ and <7i(l) = 0. That is 

9i(q) = I f{x)dx 

Jq 



Then 

3 Jiqf(q)dq 

As a second trial function we choose g%{q) = I — q s and the last trial function g^{q) = 
exp(— sx). Consider first the example given in ||, f(u) = u(l — u)(l + au), with a > 0. 
This form falls in the category given above for which an exact solution may be found. The 
transition from the regime of validity of the linear marginal stability mechanism to the regime 
of nonlinear behavior occurs at a = 2. The results obtained for this function are shown in 
Fig. [|. The dotted line labelled AW correspond to the upper bound 2^/sup(/(w)/w), and 
the dotted line labelled BN is the bound obtained from equation (3). It crosses c* = 2 at 
larger a. The solid line corresponds to the bound with the trial function g 3 with s = 7, 
the short dashed line corresponds to the bound obtained using g>2 with s = .5 and the long 
dashed line is the bound using g\ calculated form equation (7). Aronson and Weinberger's 
criterion shows that linear marginal stability is valid for < a < 1, and our bound indicates 
that it is not valid for a > 3.6. As we said above the exact solution for this case is known, 
the transition value from linear to nonlinear marginal stability occurs at a = 2. Next we 
apply the bound to the quartic polynomial / = x (1 — x)(l + ax 2 ). In this case the exact 
solution is not known and neither is the transition value from the linear to the nonlinear 
regime. The results are shown in Fig. 0, where we have used the same labelling and type 
of line as in Figure 1. Aronson and Weinberger's criterion guarantees that linear marginal 
stability is valid for < a < 4, and of the simple bounds calculated here the best shown in 
the picture corresponds to that obtained with g 2 for s = .1 which shows that linear marginal 
stability is not valid for a > 10.3. One could of course attempt to obtain a sharper estimate 
by choosing better trial functions, but this is not our purpose here. 

In conclusion, it is evident from the present results that for all non concave functions 
f(u) the linear speed cl is the asymptotic speed in a rather limited region. There is no 
substantial difference in the behavior of arbitrary polynomials, for which no exact solutions 
are known, with those already analyzed in the literature for which the exact solution and 
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point of transition can be calculated. Once the function / becomes sufficiently large, the 
selected speed will be that of the so called nonlinear front. Given the limited validity of the 
linear selection mechanism for the nonlinear diffusion equation a similar situation can be 
expected for higher order equations. Moreover, since the lower bound on the speed depends 
on the integral properties of /, it is not difficult to imagine a situation where two functions 
are identical near the origin and differ significantly near u — 1. In that case it is possible 
that the asymptotic speed for one of them be the linear value and for the other the nonlinear 
value. No local analysis of the approach to u = can then predict the transition from the 
linear to the nonlinear marginal stability regime. Finally we wish to point out that the 
analysis of monotonic fronts in phase space is useful not only in the case presented here but 
for generalized diffusion equations and in higher order equations as well. For the porous 
media equation 



with /(0) = /(l) = 0, and / > in (0,1) monotonic fronts may exist for 

> 2 SoVfvhdq 

where a(u) > must be chosen so that 

h(u) = -mu m - x a\u) > in (0,1). 
For a more general equation 

ut = (<f>(u)) xx + f(u), with 0' >0 G (0,1), 0(0) = 0, 
with the same conditions on /, monotonic fronts may exist for 



u t = (u m ) xx + f(u), 



m > 1 



c > 2 



SoVfvhdq 
Jo o d( l 



where a(u) > must be chosen so that 



h(u) 



a'{u) 4>'(u) > 0. 
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The details will be given elsewhere. It has been applied by us to obtain bounds on the 
speed of certain third order nonlinear differential equations of the type which arise in crystal 
growth problems |T| and for the dispersive Kuramoto-Sivashinsky equation ||12|| . It also 
enables one to characterize the type of functions f(u) for which the exact point of transition 
to the nonlinear regime in the the nonlinear diffusion equation (1) can be calculated without 
solving the equation explicitly [12[ . 
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FIGURES 

FIG. 1. Bounds on the speed of the monotonic fronts for the exactly solvable case 
f(u) =ii(l-ti)(l + ou). The dotted lines are the upper and lower bounds on the speed obtained 
from equations (2) and (3). The solid and dashed lines show bounds obtained from equation (6) 
with different trial functions. The bound obtained with the simple trial function exp(— 7x), shown 
with a solid line, indicates that linear marginal stability is not valid for a > 3.6. The exact value 
for the transition is a = 2. 

FIG. 2. Bounds on the speed of the monotonic fronts for the quartic polynomial 
f(u) = u(l — u)(l + av 2 ) for which the exact solution is not known. The labelling of curves 
is as in Figure 1. The bound obtained with the trial function 1 — u s with s=.l, shown with the 
short-dash line, indicates that linear marginal stability is not valid for a > 10.3. 
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